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Abstract 

This paper presents a general model for calculating the density of states and the 
Cooper pair potential in proximised superconducting bi- and trilayer films. It is valid for 
any kind of bilayer S1-S2, whatever the quality of the materials Si and S2, the quality 
of the S1-S2 interface and the layer thicknesses. The trilayer model is valid for a thin S3 
layer, whereas the other two layers have arbitrary thicknesses. Although the equations of 
the dirty limit are used, it is argued that the model stays valid in clean bi- and trilayer 
films. The typical example of superconducting tunnel junctions is used to show that 
existing models, applying to very thin or very thick layers, or to perfectly transparent 
S1-S2 interfaces, are too restrictive to apply to any bilayer. The new model is applied to 
existing junctions, with layer thicknesses intermediate between the 'thick' and the 'thin' 
approximation . 
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I INTRODUCTION 



Understanding the proximity effect in superconducting films is important for the development 
of practical devices such as superconducting tunnel junctions (STJ's). Depositing a super- 
conductor Si onto another S 2 modifies the properties of both Si and S 2 materials. If both 
superconductors are thick enough (typically thicker than 10£si(,s 2 ), with ^s 1 (s 2 ) the coherence 
length of Si (S 2 )), the extremities of the bilayer behave as bulk materials obeying the BCS the- 
ory, though not necessarily like bulk Si and bulk S 2 . The intermediate region around the Si-S 2 
interface is characterised by a relatively sharp transition between the two bulk-like regions, and 
can be pretty far from a BCS-like description. If the layers are relatively thin, any BCS-like 
behaviour can be absent from the structure. Finally, in the case where the layers are extremely 
thin, as described by McMillan [1] , each layer behaves again like a BCS superconductor. 

The physical quantities affected in a proximised bulk superconductor, are the Cooper pair 
potential A, the density of states for the Cooper pairs, P, and the density of states for the 
quasiparticles, N. As the density of states in both superconductors is modified due to the 
proximity effect, the resultant bandgap A g lies at an intermediate value between the bulk 
values for Si and S 2 , A g ,Si an d ^g,s 2 respectively. This feature has been fully described in the 
specific case, of a thin, low bandgap material S 2 next to a thick, high bandgap material Si, 
with both superconductors in the dirty limit [2, 3]. 

The goal of the present paper is to present the need for, and develop a model of, the 
proximity effect, which is not restricted to this very specific case. 

In particular section V shows that there are many situations where this special case does 
not apply, and for which the simple BCS approach does not provide a satisfactory description. 
Specifically, in the case of STJ's used as photon detectors, a more general description of the 
proximity effect is required to adequately address such issues as device performance. 

The general conditions for an extended model are presented in Section III. Such a model 
can be applied not only to STJ's, but also to any application involving bi- or trilayers of 
superconductors. The single factor which would make the model inapplicable is the roughness 
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of the S1-S2 interface. The limitations caused by this parameter are therefore also discussed. 

The extended model is described in Section II, while numerical calculations based on the 
model are presented in Section IV. The model is applied to existing STJ's in Section V, so as 
to highlight a typical application of the model. 

II The generalised proximity effect model 

The proximity effect in dirty S1-S2 bilayers was studied previously in [3], [4] -[6] for the case 
of arbitrary transparency of the Si-S 2 interface, but only the case of a thin S 2 layer and bulk 
Si layer was considered. More recently, the case of arbitrary thickness of S 2 and Si layers was 
considered in [6], but only the limiting case of the high transparency of the Si-S 2 interface 
with T Cj 5 2 = was studied. Here we generalise the theory to the case of arbitrary S 2 and Si 
layer thicknesses for the most general case of a finite critical temperature for the S 2 layer with 
T C} s 2 < T Ci 5j . We also consider the case of trilayer systems S\ — S 2 — S3 with a very thin S3 layer. 
The trilayer model remains valid for any thickness of Si and S 2 . The model is constructed for 
the case of T CjSs < T CjSl , but T C:S . 2 can have any value. 

In case where the relevant length scales (coherence lengths, mean free paths and sample 
thicknesses) are much larger than the atomic scale, the proximity effect can be described in 
terms of a quasiclassical Green's function formalism [7]- [11]. 

For an arbitrary relationship between the electronic mean free path and the coherence length 
the density of states can be derived from the solution of the Eilenberger equations [8] . If the 
vector potential is zero the Eilenberger equations have the form 
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u n a 3 - A(r)di - - (g(u n , r)) , g(u n , r) 

T 



0. (1) 



In this expression, vp is the Fermi velocity, g = d 3 g + dij is the matrix Green's function, 
A is the Cooper pair potential, r is the scattering time of the quasiparticles, and B{ are Pauli 
matrices. The energy is quantified by the Matsubara frequency, given by u) n = nT(2n + 1) (T 
is the temperature); the available energies e for the quasiparticles are related to the Matsubara 
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frequencies by the relation u n = —it. The brackets (...) denote an averaging over the Fermi 
surface. In the dirty limit, the scattering centres ensure that the Green's functions are isotropic 
all through the sample. This fact allows Eq. (1) to be simplified into a diffusion equation for 
a Green's function G(u n ,r), which is independent of the angles of the vector v F . In this case 
the function g may be represented by the form [11] 

g(u n , r) = G(u n , r) + nGi(w n , r), G 1 « G, (2) 

n being the unit vector in the direction of the momentum of the quasiparticles. 

It can be shown that G\ = v f tG-^G [9, 11]. In the case of a short mean free path, the 
quantity vpr is small, and the condition G\ <C G which is required in Equation (2) is satisfied. 
This is the condition which is realised either in the case of a high concentration of impurities 
in a bulk superconductor or dislocation scattering centres, or in the case of strong diffusive 
boundary scattering in a film. 

We consider a superconducting film Si in which the dirty limit condition l$ i < is fulfilled. 
Here is the electronic mean free path. The coherence length £ Si is related to the diffusion 



coefficient D$ i by the relation = yjDs i /2TvT Cj s i - We shall define the x axis as perpendicular 
to the film surface. The Usadel equations [9], directly derived from Eq.l can be written as 

Zs/k ( x ) + ie sin % ( x ) + A ^ (x) cos 9 Si (x) = (3) 
where the pair potential A# j is determined by the self-consistency relation 

0. (4) 



A Si (x)ln^- + 2TY: 

C,Si U) n 



A s< (rr) 

sm 6 Si {iuj n ,x) 



Here the function 6^ has been introduced as a unique Green's function defining the quasiparticle 
density of states N Si : 

N s Xe,x)/N Si (0) = Re [cos %(e, x)] (5) 

where N Si (0) is the electronic density of states in the normal state at the Fermi surface. 
In a multilayer structure the Usadel equations Eqs.3, 4 must be solved in each layer with 
the use of the appropriate boundary conditions. For convenience we substitute the coherence 



length £s t by the quantity = t;s t \jT Cy sjT C)Sl , normalised to the critical temperature of the 
layer Si. 

A Bilayer S x - S 2 

Let us first consider the case of a bilayer. For any Si-S 2 system with T C>S2 < T c ^ Sl , the origin 
of the coordinates x = is taken at the Si-S 2 interface. The region with x > is Si and 
x < corresponds to S 2 . The layers Si and S 2 have a thickness ds 1 and ds 2 , respectively. The 
boundary conditions can easily be written in terms of Os 1 (s 2 ) [4]- At the Si-S 2 interface we have 

IbnCs/s. = M0 Sl - 0s 2 ) (6) 

7& 2 ^ 2 =&A- (7) 
At the free interface of both Si and 5*2 layers, the conditions are: 

6' S2 (-d S2 )=0 ; e' Si (d Sl ) = o. (8) 

The parameters 7^ and 7 involved in the boundary conditions at the Si-S 2 interface are 
given by 

These parameters can be understood as follows: 7 is a measure of the strength of the proximity 
effect between the Si and S 2 metals, whereas 7^ describes the effect of the boundary trans- 
parency between these layers. Here ps 1 ,s 2 are normal state resistivities and Rb is the product 
of the resistance of the Si-S 2 boundary and its area. 

As is shown in Ref.[3], in the case of a thin S 2 layer where ds 2 /£,* S2 1, the parameters 
defining the proximity effect are r y m = a r yds 2 / '£5 and 7^ = ®lBds 2 /£,s 2 - Here a is a correction 
factor which is a function of the ratio of critical temperatures T C: s 2 /T C:Sl , and of the thickness 
of the layer S 2 . The dependence of the parameter a on T c ^ 2 /T CtSl can be found in Ref.[3]. 
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B Trilayer Si-S 2 - S 3 

We consider now the trilayer situation. This case is rather straightforward to solve with only 
few conditions. We first assume ds 3 <C .Again the origin of the coordinates x = is taken at 
the S1-S2 interface and the S 3 layer occupies the region — ds 2 — ds s < x < — ds 2 - As mentioned 
in [2], this approximation directly implies that pair potential and density of states in S3 are 
constant through its whole thickness. The solution in the thin S3 film is [2] 

= sin9 s (-d s )+ lBN2 A S3 
cos 9s 2 ( -ds 2 ) +7bn 2 v 

The boundary conditions at the interface Si-S 2 and at the free surface of Si can be derived 
from the general bilayer problem described above 

1b Ni Cs 2 0's 2 = sMfci - 8*) (11) 
7iOA=&A- (12) 

e' Sl (d Sl ) = o. (13) 

An additional boundary condition can be introduced at the S2-S3 interface, based on the 
fact that d Ss is very small 



Cs/ S2 (-ds 2 ) = 7 7 2o;{sinM-d Sa )-A S 3} (14) 

{l + 1 2 bn 2 ^ 2 + A| 3 ) + 2 lBN2 cose S2 (-d S2 ) [u + A 53 sine S2 (-d S2 )M} 

Here the parameters and 7^ are given by 

RBi . _ PS&. 



iBNi = —- ; 7i = e ■ (15) 

Again, the coherence lengths in S2 and S3 are normalized to T Cj 5 1 . Using a similar approach 
to those adopted previously, it is also straightforward to study the case of a very thin Si layer 
with arbitrary S2 and S3 film thicknesses. 

6 



Ill Applicability of the model to the clean limit 

The main restriction of the model presented here is the use of dirty superconductors. It is 
straigthforward to show that, in most cases, the clean limit does not need to be considered. 
The difference between clean and dirty limits lies in the presence of impurities or crystallo- 
graphic dislocations in a dirty superconductor, acting as scattering centres for quasiparticles. 
If scattering centres are on average separated by a distance smaller than the coherence length, 
then the dirty limit applies. In other words, the dirty limit localy applies to a spherical region 
of radius equal to the coherence length, and centered at a scattering centre. 

The physical presence of interfaces ensures the presence of scattering centres. As long as 
the interfaces affect several atomic layers of both superconductors, the boundary conditions are 
those of the dirty limit. Moreover, over a scale length equal to the coherence length in both 
films, any solutions for the pair potential and the densities of states are also governed by the 
same equations valid in a dirty superconductor. At depths larger than the coherence length, the 
fluctuations in pair potential and in densities of states are not significant, and the behaviour 
tends towards that of a bulk superconductor [2, 3]. Anderson's theorem [12] implies that in the 
case of a bulk superconductor, the solutions in the clean limit are the same as in the dirty limit. 
Hence, any solution for the pair potential and the density of states in regions of the films deeper 
than the coherence length, is expected to be a smooth interpolation between two dirty regions 
(in the case of a clean region separated by two rough boundaries), or a smooth extrapolation 
out of a single dirty region (in the case of a clean region located between a rough boundary and 
a flat one). The solutions related to such smooth interpolations or extrapolations are identical, 
independently of the use of the equations valid for the clean or for the dirty limit. 

The clean limit only need to be considered in cases of atomically sharp interfaces, where the 
roughness a, of the S1-S2 interface is of the order of or smaller than the interatomic distance of 
both Si and S 2 (typically around 4-5 A). This case is extremely restrictive, and also technically 
difficult to achieve, even with the most advanced methods of epitaxial thin film depositions. It 
leads to the conclusion that the clean limit does not need to be considered. 
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IV Numerical results 



In this section some typical examples of pair potential and density of states for the quasiparticles 
are determined and discussed for bi- and trilayers. 

A Bilayer Si - S 2 

For this bilayer, we have solved numerically the Usadel equations (3)-(4) by a selfconsistent 
procedure similar to that described in Ref. [3], but without the approximation ds 2 /Q 2 ^ 1- 
We start from the trial pair potentials A Sl ( S . 2 )(x) and find the solutions for Os 1 (s 2 ){ UJ m x ) m the 
Matsubara representation (e = iu> n ). Using these solutions, the new pair potentials A Sl (s 2 )(x) 
are found from the selfconsistency equation. The iterations are repeated until convergence is 
achieved. Next, after the pair potentials are determined, we solve the equations (3)-(4) on the 
real energy axis e. According to Equ.(5), this method provides spatially and energy resolved 
densities of states in both layers. 

In order to demonstrate the validity of the approach and to be representative of the devices 
discussed in Section V, the cases of an Al-Nb and an Al-Ta bilayer (T Cj ai = 1-18 K, T c ^ a = 4.5 
K, T CjiVb = 9.25 K) at T = 0.3 K are considered. 

1 Al-Nb bilayer 

The parameters for this simulation are 7 = 1.3, jbn = 2.7, d^i = 1-7£ai an d dm = 4.35 £jv&- 
This choice of parameters is justified in Section V. The results are represented in Fig.l. The 
quantities represented here are the pair potential as a function of position in the bilayer (Fig. 
la) and the density of states at different positions (Fig. lb). The density of states of Fig. lb 
is shown at the free interface of Si (Nb) and of S2 (Al) (solid lines), and on both sides of the 
Al-Nb interface (dashed lines). In Al, the density of states is peaking at an energy slightly 
higher than the gap, and then roughly decreases to 1 for infinite energies. Getting closer to the 
Nb film, a smaller peak is also visible around the energy gap of pure Nb. In Nb, the density of 
states is strongly peaking around Am , and goes to 1 at infinite energy. Below Ajv&the density 
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of states is very much depressed as compared to the density of states in Al. However, it stays 
finite and importantly the energy gap is the same as in the Al film. 

The position dependence of the energy gap is represented on Fig. la (dashed line in the 
middle) . Clearly, it is not position dependent at all, despite the very strong fluctuations of the 
pair potential. It must be stressed that, even if Fig. lb suggests a small density of states in 
Nb below A Nb , this level actually corresponds to a significant number of states. As a matter 
of fact, the density of states at the free interface of Nb (solid line in Fig. lb) at 0.3Ajvb is 
~ 1.43 • 10 6 states/meV//im 3 . Despite the fact that the gap is constant, the reduction in the 
DOS still has an effect on quasiparticle dynamics, since the quasiparticles in the energy range 
0.3 A Nb < e < A Nb entering Nb are partially Andreev reflected. As shown in Ref. [13] the 
probability to enter the reduced density of states region is proportional to N(e). 

2 Al-Ta bilayer 

For this simulation the parameters are 7 = 0.05, ^bn = 3., cIai = Cu an d dr a = 1-18 £ra- Again 
the choice of this set of parameters is explained in Section V. The results are represented in 
Fig 2. The value used for ^bn is very similar to that used for the Nb-based sample, whereas 
the value used for 7 is about 20 times lower. This has a direct impact on the shape of the pair 
potential (Fig. 2a) and the density of states (Fig. 2b). The discontinuity at the Ta/Al interface 
is extremely strong. The behaviour of Ta is almost bulk-like, with a constant pair potential 
and a steep peak of the density of states at around A^ a . However there is still a considerable 
amount of states available in Ta below Aq" (typically 0.1 x 2iV =8.1 10 6 states/meV//im 3 ). 
The fluctuations of the pair potential and the density of states within Al are also rather small, 
though a clear contribution from Ta remains visible at A^ a . Despite the very weak coupling 
between Al and Ta, the gap is still constant through the whole device thickness. Such a strong 
discontinuity at the interface emphasizes the role of quasiparticle confinement in Al, due to 
Andreev reflections. 
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B Trilayer Si-S 2 - S 3 

Next we consider two representative examples of a trilayer system: Al-Nb-NbN and Al-Nb-Ta, 
which also do have practical applications. The inclusion of a thin NbN passivation layer on the 
top electrode of an Al-Nb STJ in place of the natural niobium oxides reduces the quasiparticle 
loss rate, thereby enhancing the probability of multiple tunnel processes [14]. On the other 
hand the deposition of a thin Nb layer on top of the Al facilitates the deposition of the top Ta 
in Al-Ta STJs. The best Ta based devices, in terms of detection of optical photons, have a top 
electrode of this kind. 

In both cases the Al layer has a thickness (Lai = 0.1£^ , for which the thin layer ap- 
proximation works reasonably well, while no limitation on the thickness of other layers was 
introduced. The method of numerical solution in this case is similar to that for a bilayer with 
additional boundary conditions (14), as described in the previous section. 

1 Al-Nb- NbN trilayer 

In the case of Al-Nb-NbN we have taken thicknesses d Nh /^* Nb = 3. and d>NbN/£,NbN = 3, and 
T c ,NbN = 14K". The parameters of the Al-Nb interface are 72 = 1.3 and r )BN 2 — 2-7, as in the 
previous example. For the Nb-NbN interface 71 = 2 is expected from the resistivity ratio, while 
for the barrier transparency parameter we have assumed jbnx = 1- The results of calculations of 
the spatial distribution of the pair potential in the Al-Nb-NbN trilayer are shown in Fig. 3a. The 
dashed line shows the energy gap, which is again constant across the trilayer. The densities 
of states are shown in Fig. 3b at the free interface of NbN, x = 3£ NbN , (solid line), on both 
sides of the NbN-Nb interface, x = ±0, (dashed lines), on the Nb side of the Nb-Al interface, 
x = — 3£tv& (solid line) and in Al (dotted line). Again, in Al the density of states is peaking at 
an energy slightly higher than the gap. The magnitude of the gap is higher than in the case of 
Fig.l due to the smaller thickness of Al and the influence of the NbN layer. Note, the gap in 
Al is now as big as about half of the gap of pure NbN, Ajvbjv- The density of states at the free 
interface of NbN, x = 3£jvwv 5 peaks around A NbN , while a lot of states are also present in this 



10 



layer at e < Aatwv- 

2 Al-Nb-Ta trilayer 

In the case of Al-Nb-Ta we have taken thicknesses d Nb /^* Nb = 0.5 and dxa/^Ta = 3, and 
T c ,Ta = 4:A7K. The parameters of the Nb-Ta interface are assumed to be 71 = 0.3 and 
Ibn-l = 3.8. For the Al-Nb interface 72 = 1.3 and r )BN 2 — 2.7 were chosen again. The 
distribution of the pair potential, the energy gap and the local densities of states are shown in 
Fig.4a,b. The notations in Fig.4b are the same as in Fig. 3b with NbN and Ta interchanged. As 
is seen from comparison of Figs. 3 and 4, the energy gap in the Al-Nb-Ta trilayer is lower than 
that in the Al-Nb-NbN. This is expected, since the only difference between these two cases is 
that the high gap material, NbN, is substituted by the low gap one, Ta. Moreover, the peaks 
of the density of states in Nb are sharper in the case of Al-Nb-Ta trilayer. The reason is the 
more homogeneous distribution of the pair potential in this case, as is seen from Fig. 4a. 

V Application 

The STJ would represent a typical but not exclusive example of an application of multi-layered 
superconducting thin films. Such devices used as photon detectors have produced good results 
in terms of energy resolution, charge output and quantum efficiency, for photon energies from 
the near infrared to X-rays [15]- [20]. The most important results currently achieved can be 
summarised as follows: Ta/Al junctions provide a typical response of up to 10 5 electrons per 
eV of the detected photon, at a temperature of 300 to 400 mK; the associated measured energy 
resolution is within a few percents of the predicted theoretical limit; the quantum efficiency 
at UV wavelengths is about 60%; the longest wavelength detected is currently 2/im, with the 
possibility to extend this limit to 10/xm; the same junctions operate also at x-ray energies as 
high as 6 keV, with a quantum efficiency of about 10% and an energy resolution of about 40 eV 
FWHM. Very encouraging results have been also achieved with Nb/Al [20] and NbN/Nb [14] 
STJ's. At x-ray energies (6 keV), the best results currently achieved yield an energy resolution 
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of 12 eV for 5.9 keV photons, with an Al-based junction working at 50 mK [21]. 

A major problem to be addressed is the ability to predict the junctions response to photon 
absorption. This response can be modified by changing such variables as film thickness, film 
quality, and overall device dimensions. Especially, the performance of existing devices is still 
difficult to predict when the proximity effect plays a significant role. In particular, a significant 
difference in junction behaviour has been observed when irradiated by x-rays, as opposed to 
optical photons ([20]). Note while the trend in X-ray responsivity reflects the generally increas- 
ing role of quasiparticle self-recombination with increasing photon energies, the details of this 
response have been difficult to model using existing theories. Also, adressing adequately the 
problems of energy non-linearity ( [22]- [24]) and quasiparticle trapping ([25]) requires a detailed 
knowledge of the proximity effect theory. 

A Applicability of previous models 

The aim of this section is to show that the existing models for the proximity effect are for very 
specific cases and regimes only and need generalisation in order to explain the performance of 
an STJ. 

The type of junctions in which the energy dependence of the responsivity is difficult to 
account for are described in [26]. They essentially consist of symmetrical Nb/Al/AlO^/Al/Nb 
junctions deposited on super-polished i?-plane sapphire, with 100 nm of Nb and 15 to 120 nm 
of Al. The base electrode is epitaxial, whilst the top is polycrystalline. The polycrystalline film 
is usually covered by a thin NbO layer, followed by a Nb 2 5 layer (~ 5 nm in total), due to 
the normal exposure to air. Note, the bandgap is much smaller in Al than in Nb. The Al layer 
represents the S2 layer, and the Nb the Si layer. 

A typical value for the coherence length in bulk, clean Al at K is £ai(0) ~ 1.6 /im. For the 
polycrystalline film (superscripts p indicates polycrystalline), the electron mean free path l p Al is 
limited by the grain size. Using a Transmission Electron Microscope (TEM), the grain size in 
the polycrystalline Al film has been estimated at about 40 nm. In the dirty limit, the coherence 
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length is given by Cm = ^AiU{0)/3 ~ 147 nm, and Cm = Cai^ T C)Al / T c , Nb ~ 53 nm. Thus 
the condition of the dirty limit for the polycrystalline films Si and S2 is not actually valid. 
Furthermore, this coherence length is precisely in the middle of the range of Al thicknesses 
available, and the condition of small Al thickness is also not satisfied. As for the epitaxial 
film (superscripts e), the mean free path l e Al is not limited by the grain size anymore, but it is 
simply constrained by the boundary scattering due to the thickness of the electrode (220 nm) 
and thus is always far lower than (,ai(®)- This conclusion is consistent with RRR (resistance 
residual ratio) measurements, which lead to an average mean free path of ~ 150 nm over the 
whole electrode. The coherence length Cai can therefore be assumed to be governed by a dirty 
environment, and will not exceed Cm \J^aiI^ai < 120 nm. This is not very large compared to 
either l e Al or d A \- As a consequence, in the epitaxial film, none of the assumptions associated 
with the dirty limit and with a thin Al layer are valid. 

In Nb, £jvb(0) ~ 38 nm. This is very close to the grain size in the polycrystalline film, 
implying that the Nb polycrystalline film is in between the dirty and the clean limit. If the 
layer S\ was assumed to be dirty, the coherence length is ~ 23 nm. In reality, it is somewhere 
between 23 and 38 nm. In any case, the dirty limit does not apply. In the epitaxial layer, RRR 
measurements have determined an electron mean free path l e Nb ~ 150 nm; here we are clearly 
in the clean limit, constrained by boundary scattering, such that for a Nb film of thickness 
100 nm (applicable for all examples described herein), the condition of a thick Si layer is not 
fulfilled. 

B Application of the generalised proximity effect model 

The model has been applied to two different STJ layups. The first one is a symetrical Nb based 
STJ with lOOnm of Nb and 90 nm of Al. The second is a symetrical Ta based STJ with lOOnm 
of Ta and 55 nm of Al. In order to determine the density of states and the order parameters, all 
the input parameters for the proximity effect model have to be established. These parameters 
are the critical temperature, the ratio between the thickness of each layer and the coherence 
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length, and the interface parameters. 

The coherence length in a non bulk configuration can be expressed in terms of the bulk 
coherence length £ and of the electron mean free path /. If the superconductor is in the clean 
limit, this relation is 

111 

r^ + r (16) 

In the dirty limit, where / << £, one has 



Assuming the clean limit in epitaxial base films and the dirty limit in polycrystalline top 
films, limiting the mean free path to the grain size (40 nm) in polycrystalline films, the values 
of £/vb ~23 nm, £ai 147 nm and £r a ~85 nm were found. These values are actually averages 
between the values found for polycrystalline and epitaxial films. Using these values we determine 
the ratios of thickness to normalised coherence length as used for the numerical simulations in 
Sections IV A 1 and 2. 

The method chosen to determine the values of the interface parameters is to compare an 
experimentally determined value with the theoretical estimate, as derived from a simple analyt- 
ical expression for this value, which includes the densities of states. The quantities to be fitted 
to the model were chosen to be the energy gap as a function of temperature and the critical 
current as a function of temperature. Both quantities are easily found experimentally from the 
junctions I-V curve characteristics. The critical current I c through the junction is calculated 
in the following way [3] 

TTT^ = ^fY. sin^(w n )sin^K), (18) 

where R n is the normal-state resistance of the junction, and 9i(2)(w n ) is the local Green's 
function in the vicinity of the barrier in electrode 1 and 2 respectively. Both electrodes are 
assumed to be identical, therefore simplifying Eq. (18), since Fi — F 2 , and thereby ensuring a 
unique set of (7, 7bjv) parameters for the whole junction. 
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Figs. 5a and 6a show the model values compared to the experimental data of the energy 
gap as a function of device temperature for the Nb and the Ta based STJs respectively. Figs 
5b and 6b show essentially the same information, but this time in regard to the critical current. 

The fits to the experimental points correspond to model values with 7 = 1.3, ^bn = 2.7 for 
the Nb based STJ and 7 = 0.05, ^bn = 3 for the Ta based STJ. In both cases the fit of the 
bandgap is satisfactory, while the critical current of the Nb based device is slightly too low. 

Using these model derived values of the parameters we obtain the results presented in 
sections IV A 1 and 2. 

C Dependence of the interface parameters on film thickness 

We will now analyse how the interface parameters vary when the thicknesses of the films change, 
all other parameters staying the same. Assuming that the dirty limit is valid, we can replace 
the coherence length by its expression in the dirty limit: 

t = )[f (19) 
This gives the following expressions for the interface parameters: 



1BN = C, 



1BN 1 



where the constants C 7 and C 7SAr are defined by: 



— ■ ■» = C ^M ■ < 2 °> 



r — ^ B / _ ■ r — eih. H°Z ro-n 

These two constants are independent of the film thickness. C 7 depends only on the nature 
of the two films involved, whereas C 7BJV depends also on the macroscopic properties of the 
interface. The only variables depending on the film thickness are the critical temperature and 
the mean free path. 

The critical temperature as a function of film thickness can be determined according to the 
model of Cooper [28] . This model states that superconductivity is lost in a thin surface layer b t 
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due to a reduction in the electron density of states near the surface. The critical temperature 
dependence according to this model is: 

r -- r <4-:&)- (22) 

where T Cj0 is the critical temperature of the bulk material, N is the electron density of states 
at the Fermi level, v the bulk interaction potential and t the film thickness. 

The bulk mean free path at low temperature 1 (just above the critical temperature) can be 
found via the bulk mean free path at 300K, lo,3oo, and the residual resistance ratio (RRRb) of 
a thick film: 

la = RRRblo,3oo- (23) 
We can then calculate the mean free path in a thin film (t << 1 ), using [29]: 

l = ^^n l j + 0.423^). (24) 

Using the interface parameters found for a Nb based STJ with 45 nm of Al and the basic 
material constants for Nb and Al (table 1) appearing in equations (19)-(24), we can determine 
the interface constants C 7 and C lBN . Except for the factor Nz/ (cfr. eq. 22), which was 
determined in Ref. [30], all values in table 1 were derived from experiments on Nb thin films 
or Nb/Al/AlOx/Al/Nb multilayers. For these films the critical temperature and RRR were 
measured as a function of film thicknesses. A fit to these experimental points provided the 
values of b t and T C)0 . 

Using equations (22) and (24) we then calculate the values of the interface parameters 
for different Al film thicknesses. The theoretical and experimental values of the interface 
parameters for symmetrical Nb/Al STJs with Al film thicknesses of 15, 30, 45 and 90 nm and 
a constant Nb film thickness of lOOnm can be found in table 2. The experimental values were 
determined in exactly the same way as described in the previous section. The correspondence 
between experiment and theory is very good, when we take into account the rather large 
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uncertainties we are dealing with respect to the numerous parameters involved. Note that 
including the •jbn parameter in the discussion is important, as it is not small and has a strong 
effect on the density of states in the two films. This is in accordance with the results obtained 
by Zehnder et al. ( [31]). Nevertheless the dependence on film thickness is rather complex for 
both interface parameters. No clear square root dependence for 7 or linear dependence for jbn 
was observed, as stated in ref [31]. 

VI CONCLUSIONS 

A model describing the proximity effect for any kind of superconducting bilayer, in terms of 
thickness, critical temperature and cleanliness, and superconducting trilayers with a thin third 
layer has been presented. It has been proposed that any existing samples, even if very clean, 
would obey to these equations which are valid in the dirty limit, because of the presence of 
imperfections at the boundary. Only in very specific cases of extremely thin and smoothly de- 
posited layers (with an rms roughness at the interface of the order of the interatomic distance 
in both materials) would the model not apply. This model has been presented and examined 
experimentally for typical values of the parameters involved for tunnel junctions, using thick- 
nesses intermediate between the extreme cases discussed in previous publications. Finally, the 
model has been shown to be very effective in determining the various important parameters 
for practical cases, using the current-voltage characteristics of STJ's (energy gap and critical 
current as a function of temperature). 
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VII Tables 
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Table 1: Basic parameters for Nb and Al 





7 (th.) 


7 (exp.) 


7B7V (th.) 


7btv (exp.) 


15nm Al 


1.07 


0.9 


1.85 


1.5 


30nm Al 


1.32 


1.2 


2.3 


2.2 


45nm Al 


1.6 


1.6 


2.7 


2.7 


90nm Al 


1.6 


1.3 


2.7 


2.7 



Table 2: Comparison of theoretical and experimental values of the interface parameters for 
Nb/Al STJs with different Al film thicknesses. 
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Captions for figures 
Figure 1: 

Typical output for Al-Nb. Pair potential (a) and density of states (b) for 7 = 1.3 and 
•Jen = 2.7, with ds 2 = 1-7 £2 an d ds l = 4.35 £1. (a) The upper dashed line is the bulk energy 
gap of Nb. The lower dashed line is the bulk energy gap of Al. The intermediate dashed 
line is the resulting energy gap, as determined from Fig 1 (b). (b) The densities of states 
are represented for both materials at the free interface (solid lines) and at the S1-S2 interface 
(dashed lines). 

Figure 2: 

Typical output for Al-Ta. Pair potential (a) and density of states (b) for 7 = 0.05 and 
Ibn — 3., with ds 2 = £2 an d ds 1 = 1.18 £1. (a) The upper dashed line is the bulk energy gap 
of Ta. The lower dashed line is the bulk energy gap of Al. The intermediate dashed line is the 
resulting energy gap, as determined from Fig 2 (b). (b) The densities of states are represented 
for both materials at the free interface (solid lines) and at the S1-S2 interface (dashed lines). 

Figure 3: 

Typical output for Al-Nb-NbN. Pair potential (a) and density of states (b) for 71 = 2., 
Ibn! = 72 = 1-3 and r yBN 2 — 2.7 with ds 3 = 0.1 £3, ds 2 — 3 £2 an d ds 1 = 3 £1. (a) The upper 
dashed line is the bulk energy gap of NbN. The lower dashed line is the bulk energy gap of Al. 
The intermediate dashed line is the resulting energy gap, as determined from Fig 3 (b). (b) 
The densities of states are represented for NbN at the free interface (solid line), at the S1-S2 
interface (dashed lines)and at the S 2 — S 3 interface (dotted lines). 

Figure 4: 

Typical output for Al-Nb- Ta. Pair potential (a) and density of states (b) for 71 = 0.3, 
7 BJVl = 3.8, 72 = 1.3 and 7bjV2 = 2.7 with d Ss = 0.1 £*, d S2 = 0.5 £ 2 * and d Sl = 3 £ x . (a) The 
upper dashed line is the bulk energy gap of Nb. The lower dashed line is the bulk energy gap 
of Al. The intermediate dashed line is the resulting energy gap, as determined from Fig 4 (b). 
(b) The densities of states are represented for Ta at the free interface (solid line), at the S1-S2 
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interface (dashed lines) and at the S 2 — S 3 interface (dotted lines). 
Figure 5: 

Energy gap (a) and critical current (b) as a function of temperature for the Al-Nb STJ 
(cIai = 90nm, dm = lOOram). The data (diamonds) is fitted using 7 = 1.3 and •jbn = 2.7.The 
expected experimental uncertainty for the data points is also represented. 

Figure 6: 

Energy gap (a) and critical current (b) as a function of temperature for the Al-Ta STJ 
(d,Ai = 55nm, dr a = lOOnm). The data (diamonds) is fitted using 7 = 0.05 and jbn = 3. 
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Pair potential for Al-Nb 
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